Morphogenesis of membrane invaginations in spherical confinement 
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We study the morphology of a fluid membrane in spherical confinement. When the area of the 
membrane is slightly larger than the area of the outer container, a single axisymmetric invagination 
is observed. For higher area, self-contact occurs: the invagination breaks symmetry and deforms into 
an ellipsoid-like shape connected to its outer part via a small slit. For even higher areas, a second 
invagination forms inside the original invagination. The folding patterns observed could constitute 
basic building blocks in the morphogenesis of biological tissues and organelles. 

PACS numbers: 87.16.D-, 87.10.Pq 



Introduction. Folding phenomena are ubiquitous in 
nature and living matter. They are a key to understand- 
ing the complex shapes of the mammal gut [1 , the cere- 
bellum [2] and the kidney [3 j, and are of fundamental 
importance in developmental biology. 

Typically, folding processes are triggered by either a 
reorganization of mass due to growth [4] |5j [8] or by a 
buckling instability originating from external forces or 
constraints [6j[7][9]- Here, we study the folding deforma- 
tions of a thin closed membrane inside a spherical cavity. 
This system can be considered a prototype for surface 
invagination, a process that occurs in different biological 
systems. Among others, well-known examples are the 
gastrulation within the surrounding egg shell [10] or the 
crista formation of the inner mitochondrial membrane 
inside the outer membrane [IT]. In both cases, the fold 
formation is a geometric necessity, since the inner surface 
is too large to fit into the cavity. Scientists have studied 
the resulting shapes and how they develop during growth 
since decades. 

In this Letter, we use a minimal mechanical model of a 
fluid membrane to study the shape of such invaginations. 
Using finite element simulations, we construct a morpho- 
logical phase diagram showing what shapes emerge for 
given membrane area and enclosed volume. We find an 
initial invagination of axisymmetric type for small surface 
growth, which breaks symmetry into ellipsoidal shapes as 
the membrane surface is increased. For even higher pre- 
scribed area, a secondary invagination emerges inside the 
first one as a result of the interplay of volume and sur- 
face constraints together with self-contact. Our analysis 
is capable of showing the morphogenesis of cristae-like 
invaginations with as few ingredients as possible, and ex- 
tends previous works of, e.g., Ref. [12 considerably. 

The model. The classical curvature model developed 
by Canham, Helfrich, and Evans expresses the mechan- 
ical energy of a membrane related to its bending as a 
second order expansion in curvatures [T3HT5] : 

E b =^ I dA(2H-C ) 2 + K G f dAK , (1) 



where the integrals are carried out over the surface of 
the membrane. Here, H is the mean curvature and K is 
the Gaussian curvature. The spontaneous curvature Co 
represents an intrinsic preferred mean curvature of the 
membrane which we will set to zero for simplicity, k and 
kq denote bending and Gaussian rigidity, respectively. 
For a given topology, the second term in the energy ([I]) 
equals a constant due to the Gauss-Bonnet theorem. We 
can omit this term here since only closed membranes with 
spherical topology are considered. Fixing area and vol- 
ume to A and V, respectively, the scaled total energy of 
the membrane can thus be written as 

E=[dA2H 2 + ^(A-A) 2 + ^(V-V) 2 , (2) 
J z z 

where \±a and jly are penalty factors implementing the 
constraints. This model is not only relevant for fluid lipid 
bilayer membranes but can also be applied to growing 
soft tissue as long as in-plane shear forces can equilibrate 
on time scales smaller than growth [16]. 

We mimic the effect of the external constraint by 
adding a soft spherical container to the system. This con- 
tainer is modeled via a spherical force field which penal- 
izes the motion of the membrane outside the container. 
Applying such a soft constraint implies that the mem- 
brane is allowed to trespass into the force field at the 
expense of increasing its energy. As a result, the radius 
of the overall shape can be slightly larger than the radius 
of the outer shell during the simulations. 

Finite element simulations. The bending energy in- 
volves squares of curvatures, which are second derivatives 
of the surface vector function. For a finite element sim- 
ulation it is possible to use either a mixed method as 
suggested in Ref. [17 or trial functions that have square 
integrable derivatives up to second order. An elegant 
formulation of such elements for thin shells has recently 
been developed by Cirak et al. The method is based on 
a CAD subdivison scheme which satisfies this C 1 conti- 
nuity requirement of the surface vector function [18] [19] . 
Klug and coworkers extended this approach to fluid lipid 
membranes [20] [21] . Following these works we discretize 
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FIG. 1: (a) Numerical equilibrium solution for a membrane 
(smooth surface) inside a spherical container (black mesh) 
with scaled area a — 1.2 and scaled volume v = 0.8. The 
membrane bulges inward and forms an invagination reminis- 
cent of a light bulb, (b) Vertical slice of the system, which 
contains the axis of symmetry, together with the correspond- 
ing bending energy density of the membrane (= 2H 2 ). In (b) 
the measured distances of the invagination are defined (see 
Tab. 



the scaled membrane energy given in eq. Q and set up 
corresponding forces that act on the nodes of the mesh. 
To prevent mesh distortions and possible numerical insta- 
bilities, we applied viscous regularization and r-adaptive 
remeshing schemes as described in [21]. Instead of using 
line search based solvers as in Ref. [2T] , we performed a 
time integration and added a damping force on the nodes 
to equilibrate the system analogous to the method used 
in Ref. [22] for thin shells. 

Moderate surface growth We simulated the membrane 
in a spherical container of unit size for different values of 
scaled surface area a = A/Aq and scaled volume v = 
V/Vo, where Aq and Vq are the area and the volume of 
the container. 

If the area of the membrane is larger than the area of 
the container, the membrane has to fold into the inte- 
rior. For moderate values of surface growth, the mem- 
brane forms a single invagination, connected via a neck 
to the part in contact with the outer shell (see Fig. [I]). 
Every additional invagination would contribute an extra 
energy of about 87r^ — the value of the bending energy 
of a sphere — to the overall bending energy which is why 
multiple invaginations are not observed in equilibrium. 

The shapes that we obtain for moderate surface growth 
appear to be axisymmetric. To confirm this, we took hor- 
izontal slices of the surface and analyzed their curvature: 
first, we determined the symmetry axis by taking the 
mean of the surface normals of each vertex. Second, we 
took a slice of the membrane perpendicular to the axis 
of symmetry at the point where the invagination is the 
thickest. Using a discrete formulation of the curvature 
based on a second order polynomial approximation of 
the curve on sample points [23], we estimated the curva- 
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TABLE I: Measurements of characteristic system parameters 
see Fig. [1(b)] . The error of measurement is about ±0.05 
for all parameters including the scaled bending energy eb := 
E b /(8tt«). 



ture for each invagination slice. The results confirm that 
the obtained shapes are axisymmetric for moderate sur- 
face growth. In this case it is thus sufficient to consider 
a two-dimensional vertical slice of the shape which con- 
tains the axis of symmetry. The corresponding slice for 
(a, v) = (1.2, 0.8) is shown in Fig. 1(b) | together with the 
bending energy density. Although the membrane shape 
resembles a stomatocyte as is found in the corresponding 
reduced volume problem without container, we note that 
it is the confinement which forces the membrane into this 
form; a free vesicle would adopt the form of an ellipsoid 
for these parameter values [24] . 

To analyze the shapes quantitatively, we examine the 



corresponding system parameters defined in Fig. 1(b) 
the variable L\ denotes the horizontal distance between 
the uppermost points, and L 2 and L 3 are taken at the 
narrowest and the broadest point of the invagination, re- 
spectively. The vertical distance h\ indicates how much 
the membrane detaches from the container; h 2 measures 
the length of the neck between the uppermost points and 
the points where the curvature of the slice changes its 
sign. The distance hs is the extension of the tip of the in- 
vagination. In Tab. [I] the corresponding numerical values 
are listed for v = 0.8 and 0.9. We remark that for con- 
stant volume v and increasing surface area a, the length 
hs of the tip increases, whereas the neck decreases in 
length h 2 and diameter L 2 . Increasing the inner volume 
from 0.8 to 0.9 does not significantly change the overall 
shape, but causes the invaginations to penetrate less (see 
Fig.§. 

We have also measured the bending energy e^ of the 
whole membrane shape, normalized by the bending en- 
ergy of a sphere (see again Tab.[l|). We find that e^ is very 
close to two, i.e., the energy of a system of two isolated 
spheres. We therefore consider such a system as a simple 
model for the membrane. In this case, there is no connec- 
tion between the invagination (i.e., the inner sphere) and 
the part of the membrane in contact with the container 
(the outer sphere). The radius of the outer sphere is one, 
the radius of the inner sphere is set to Ri. The ratio a of 
total area to the area of the container can be written as 
a = 1 + B% . Similarly, the ratio v of the volume between 
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FIG. 2: Morphological phase diagram with vertical slices of 
each simulation point. The color of each slice corresponds to 
the different regimes observed: the invagination is either ax- 
isymmetric (black), ellipsoid- like (dark gray), or stomatocyte- 
like (light grey). The slices of the non-axisymmetric shapes 
lie in the symmetry plane perpendicular to the slit-like neck. 
For the dotted slices a metastable ellipsoid-like state has been 
observed. Eq. (|3| is plotted for different values of v (dot- 
ted lines). Inset: mean curvature H of the vertical slice 
(a, v) = (1.2, 0.8) as a function of scaled arc-length s/(2n). 

the two spheres to the volume of the container is given by 
v = 1 — Rf . Hence, we get the following simple relation 
between the two parameters: v = 1 — (a — l) 3 / 2 . If this 
equation is fullfiled, both constraints can be accomodated 
by the 2-sphere system. At a = 1 and v = 1, the inner ra- 
dius vanishes and we are just left with one unit sphere. At 
a = 2 and v = 0, the inner sphere reaches the size of the 
outer one, and the system consists of two unit spheres. If 
v < 1 — (a — l) 3 / 2 , the system will be stressed, since the 
volume between the two spheres wants to be smaller than 
the membrane can accomodate with the available area. 
Naively, one would think that no axisymmetric ground 
state exists any more. In the simulations, however, the 
outer part can detach from the container. This allows 
the membrane to decrease the enclosed volume such that 
an axisymmetric state is still energetically favored. 

A visual inspection of the shapes obtained suggests 
that the neck can be approximated theoretically by a 
catenoidal batch connecting the two spheres. However, 
such a model turns out to be insufficient to explain the 
results of the simulations quantitatively. A quick glance 
at the mean curvature along a vertical simulation slice 
explains why [see inset of Fig. [2]: the mean curvature H 
at the neck is not close to zero but changes abruptly from 




(c) (d) 



FIG. 3: Numerical equilibrium solutions which break axisym- 
metry. (a) Ellipsoid-like state for (a,v) — (1.6,0.7). (b) 
Stomatocyte-like state for (a, v) = (1.6,0.9). (c) Cut of the 
stomatocyte-like state along the symmetry plane, (d) Cut of 
the stomatocyte-like invagination perpendicular to the sym- 
metry plane. 

negative (part in contact with the container) to positive 
(invagination). 

For v > 1 — (a — 1) 3 ^ 2 , too much volume wants to be 
stuffed into a membrane of limited area. One may antic- 
ipate that this frustration is resolved by breaking sym- 
metry. This happens indeed as one can see in Fig. |3(a)| 
However, the neck mitigates this effect: the tip of the 
invagination can still keep a shape close to a sphere if v 
is not too large (see Fig. [2J. 

High surface growth and symmetry breaking Packing 
increasing membrane area into a constrained volume in- 
evitably produces self-intersections. To avoid this, we 
implemented a collision handling algorithm based on the 
minimization of the intersection contour of polygonal el- 
ements [29] . 

With this algorithm, non-axisymmetric systems can be 
studied as well. For (a, v) = (1.5,0.7), and (1.6,0.6) — 
i.e., systems with an axisymmetric global minimum — a 
metastable symmetry breaking state with self-contacts 
has been observed in the simulations. At the connect- 
ing neck, the membrane touches itself and forms a non- 
spherical slit, which endures for long simulation times, 
indicating that this configuration is a local minimum. 
The invagination is not spherical in this case. Its form 
is reminiscent of the oblate and prolate shapes found in 
reduced volume problems such as in [24] . 
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To include this behaviour in the theoretical model, we 
again neglect the neck and assume that the system con- 
sists of an outer sphere of radius one (i.e., the part of the 
membrane in contact with the container) and a vesicle 
of volume uAttR^ /3 (the invagination), where the param- 
eter v G {0, 1} is the reduced volume which measures 
the deviation from the volume of a sphere. The relation 
between a and v now reads: 
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According to Ref. [24] the theoretical model adopts the 
following shapes: for 0.65 < v < 1 the global energy 
minimum corresponds to a spherical volume enclosing a 
prolate ellipsoid in the middle. Decreasing the value of 
v more and more, causes the ellipsoid to become oblate 
for 0.59 < v < 0.65, while for v < 0.59 the invagination 
distorts into a stomatocyte shape. In our simulations 
corresponding shape transitions can indeed be found (see 
dotted lines in Fig. [2]). The connection via the neck, how- 
ever, alleviates these transitions. Moreover, the confine- 
ment gives rise to more complex invagination geometries: 
we distinguish two different types of non-axisymmetric 
invaginations in the simulations. Ellipsoidal invagina- 
tions have a complex slit-like neck: the middle of the slit 
touches itself so that it appears as if there are only two 
holes as openings [see Fig. 3(a) . Increasing the surface 



area of the membrane or the volume enclosed by it de- 
creases v\ a second invagination bulges into the first one 
[see Fig. |3(b)| - 3(d)| . The shape of this secondary invagi- 
nation resembles a stomatocyte [see Fig. |3(d)| white line] 
connected to the outer part via a slit-like neck. 

Conclusions Combining a simple mechanical model 
of a closed fluid membrane with the constraint of a con- 
fining shell, we investigated the morphogenesis of mem- 
brane invaginations. A theoretical approximative model 
allowed us to classify the shapes in a morphological phase 
diagram. Owing to the constraint, the equilibrium shapes 
deviate considerably from the reduced volume solutions 
found for membrane vesicles without confinement [24] . 
The axisymmetric invaginations obtained for moderate 
surface growth resemble closely the basic geometry of bi- 
ological invaginations such as inner mitochondrial cristae 
or embryonal gastrula. Incorporating the treatment of 
self-contacts of the membrane into the model, we found 
non-axisymmetric shapes for high surface growth, namely 
ellipsoid-like and stomatocyte-like invaginations. Our 
computational analysis of a minimal membrane model 
could help to establish which structures are intrinsic to a 
folded membrane and which may be particular to other 
effects such as membrane protein interactions in the case 
of mitochondria [25] [26] or localized cell differentiation 
in gastrulation [27] . 
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